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Abstract. We have re-analyzed the 6-12/jm ISO spectrum of the ultra-luminous infrared galaxy Arp220 with the conclusion 
that it is not consistent with that of a scaled up version of a typical starburst. Instead, both template fitting with spectra of the 
galaxies NGC4418 and M83 and with dust models suggest that it is best represented by combinations of a typical starburst 
component, exhibiting PAH emission features, and a heavily absorbed dust continuum which contributes ~40% of the 6- 
12/jm flux and likely dominates the luminosity. Of particular significance relative to previous studies of Arp220 is the fact 
that the emission feature at 1.1 pm comprises both PAH emission and a broader component resulting from ice and silicate 
absorption against a heavily absorbed continuum. Extinction to the PAH emitting source, however, appears to be relatively 
low. We tentatively associate the PAH emitting and heavily dust/ice absorbed components with the diffuse emission region and 
the two compact nuclei respectively identified by Soifer et al. 1 2002 ( in their higher spatial resolution lOyum study. Both the 
similarity of the absorbed continuum with that of the embedded Galactic protostars and results of the dust models imply that 
the embedded source(s) in Arp 220 could be powered by, albeit extremely dense, starburst activity. Due to the high extinction, 
it is not possible with the available data to exclude that AGN(s) also contribute some or all of the observed luminosity. In this 
case, however, the upper limit measured for its hard X-ray emission would require Arp 220 to be the most highly obscured AGN 
known. 
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1. Introduction 

■ The galaxy Arp 220 (IC4553; cz=5450km/s; D=73Mpc) was 
originally classified by Arp J1966i as a "galaxy with adjacent 
loops". Its optical image (1" = 352 pc) shows faint structures, 
reminiscent of tails or loops, suggesting it to be the remnant of 
a recent galaxy merger (Toomre & Toomre [l972i ). IRAS (1983) 
increased interest in Arp 220 through the discovery of its far-IR 
luminosity and infrared-to-blue ratio which characterized it as 
an extreme member of the "unidentified infrared sources" dis- 
covered during the mission (Houck et al. 119841 Soifer et al. 
I1984> . When later the spectroscopic redshifts of these "uniden- 
tified infrared sources" became available, Arp 220 turned out 
to have a similarly large infrared luminosity (1.35 x 10 12 L G ), 
making it the nearest member (by a factor of ~2) of the new 
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class of UltraLuminous InfraRed Galaxies (ULIRGs; Sanders 
et al. 119881 . with L iR >10 12 L Q . Numerous studies across all 
wavebands have since examined Arp 220 in close detail, also 
showing this nearest ULIRG to be unusual in some aspects 
rather than being typical for the class. 

Like most other ULIRGs, Arp 220 is the product of the 
interaction of two gas-rich disk galaxies (Sanders & Mirabel 
1996). Groundbased observations at 10-30yum suggest that its 
luminosity arises in the innermost 250 pc (Wynn- Williams & 
Becklin 1993). Radio and mm observations reveal its two nu- 
clei to be surrounded by molecular disks of r~100pc, which 
counterrotate with respect to each other (Sakamoto et al. 1999 ). 
The eastern nucleus seems to be embedded within an outer gas 
disk of r~l kpc, which rotates in the same sense. The western 
nucleus is connected to the eastern nucleus by a thin gas bridge, 
traced in H I absorption, and appears to lie above the outer gas 
disk (Mundell et al. 2001 ). The projected separation of the two 
nuclei amounts to 345 pc (0.98"; Baan & Haschick [T995l ). High 
sensitivity VLBI observations disclose the presence of multi- 
ple compact radio sources dispersed over the two nuclei. The 
knots are consistent with free-free emission from luminous ra- 
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Fig. 1. The 2-3000 fim spectrum of Arp 220. Filled circles represent the IRAS fluxes. Spectra shown are: the smoothed 2.4-4.9 fim 
ISO-PHT-S spectrum (Spoon et ai. l2"002l : the 3.2-3.8/mi CGS4 spectrum obtained in 1.2" slit (Imanishi & Dudley l2000b : the 



5.0-16/mi ISO-CAM-CVF spectrum (Tran et al. 200D; the 17-22/mi UCL spectrum (Smith et al. 1989); the smoothed 45- 
200 fim ISO-LWS spectrum (Fischer et al. 1997 1. Other spectral data included in the plot are ISO-PHT photometry (Klaas et al. 
119971 M.Haas, priv. comm.); ISO-SWS background subtracted continuum measurements (E.Sturm, priv. comm.); UKIRT and 
SCUBA far-IR photometry (Eales et al. [T989l Dunne et al. l2000ll200Tl and mm-observations (Anantharamaiah et al.i20()0J. The 
dash-dotted line is one choice for the local continuum in the 5-25 /mi region (see also Fig.[5}. The inset shows a comparison 
of the 5.8-1 1.7/rni ISO-PHT-S (Spoon et al. 12003 and the 5.0-16/mi ISO-CAM-CVF spectra (shown as grey and black lines 
respectively) with the Keck-MIRLIN photometry (filled circles) of Soifer et al. J1999I I. 



dio supernovae expanding in a dense medium (Smith et al. 
119951 119991 . 

At shorter wavelengths (in the UV, optical and near-IR) 
the view towards the nuclear components is greatly impaired 
by strong dust extinction of at least A v =30-50 (Sturm et al. 
119961 . In the mid-IR, the dust opacity (A A ) is a factor of 10-100 
less than at optical wavelengths and Smith et al. 1^989) used 
this property to study the nature of the central power source 
in Arp 220 in the 8-1 3 /mi (N-band) and 17-22 /mi (Q-band) 
atmospheric windows. Based on the weakness of the 11. 2 /mi 
PAH emission band within the deep 9.7 fim silicate absorption 
feature, they concluded that only 2-10% of the total infrared 
luminosity is powered by starburst activity, with an obscured 
AGN responsible for the rest. Further analysis of the 9.7 fim sil- 
icate absorption feature led Dudley & Wynn- Williams ( 1997) 
to conclude, however, that the obscured power source resem- 
bles a scaled-up embedded protostar. 

Not limited to the mid-IR atmospheric windows, ISO spec- 
troscopy revealed two pronounced spectral features in the pre- 
viously unstudied 5-8 fim range. In line with ISO observations 
of Galactic star forming regions, the two features were identi- 
fied as the 6.2/mi and 7.7/zm PAH emission bands (Genzel et 
al 1998 1. Using the ratio of 7.7 fim PAH emission to the un- 



derlying 7.7 fim continuum as a criterium to discern starburst- 
and AGN-dominated galaxies, Genzel et al. ( 1998 1, Lutz et al. 
JT998t . Spoon et al. JT998t . Rigopoulou et al. JT9991 and Tran 
et al. (2001) classified Arp 220 as starburst-dominated. High 
angular resolution groundbased N-band spectroscopy has since 
shown the 11.2 fim PAH emission in the nuclear region to be 
diffusely distributed over the central ~2" and the starburst as- 
sociated with the PAH emission not to be able to account for 
more than 10-50% of the bolometric luminosity (Soifer et al. 



2002). In summary, the mid-IR low -resolution spectral diag- 
nostics appear mostly starburst-like but star formation appears 
quantitatively insufficient to account for the bolometric lumi- 
nosity, unless strongly obscured or otherwise modified. The 
same is, to a lesser degree, true for the more direct tracing of 
starburst activity by mid-IR fine-structure lines. The ratio of 
~l/37 of starburst ionizing luminosity and bolometric luminos- 
ity, derived by Genzel et al. 



; 1998 1 for this source, is about a 
factor 2 less than in comparison starbursts. 

Since Arp 220 is often regarded as a nearby template for 
dusty galaxies at high redshift undergoing vigorous star for- 
mation (e.g. faint SCUBA sources), it is imperative to clearly 
identify its power source(s). Despite the quantitative problems 
with the starburst origin for the luminosity, alluded to above, 
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the general consensus since ISO has been massive young stars 
(Genzel & Cesarsky 2000). However, the infrared luminous 
galaxy NGC4945, which shows no outward evidence for an 
active galactic nucleus even in ISO observations (Genzel et al. 
1998 Spoon et al. 2000), has turned out to contain a heav- 
ily obscured AGN visible only in hard X-rays (Iwasawa et al. 
[E>93l Done et al. [19961 Guainazzi et al. 2000). For Arp 220, 
BeppoSAX and Chandra observations do not detect a similar 
hard X-ray source (Iwasawa et al. 2001 , Clements et al. 2002 1. 
The only possibility for an energetically significant AGN to 
exist in Arp 220 would therefore be in the form of a deeply 
embedded source, hidden behind a 'Compton-thick' shell of 
Nh >10 25 crrT 2 with a covering factor close to unity (Iwasawa 
et al. 2001 1. The presence of huge amounts of molecular gas in 
the central parts (~10 10 M o ; Scoville et al l 19971 Sakamoto et 
al. 1999 ) indicates that sufficient obscuring material is indeed at 
hand. And the very large 850 fim dust-continuum flux to 7.7 fim 
PAH flux (Haas et al. 2001 ) could mean that the luminosity of 
this embedded source is redistributed into the far-IR. 
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2. The infrared spectrum of Arp 220 

The 2-3000 fim spectral energy distribution of Arp 220 (Fig.^ 
is characterised by a steeply rising dust continuum, setting in 
around 4-5 fim, and leading up to a remarkably strong far-IR 
flux peak at 60-100 fim. The most prominent features in the 
infrared spectrum are the silicate absorption features at 9.7 fim 
and 18//m and a broad flux peak at 7.8 fim. 

2.1. The mid-IR spectrum of Arp 220 

In Fig.FJlwe compare the mid-IR spectra of Arp 220, NGC 4418 
and the central region of the starburst galaxy M 83. At first sight 
the three spectra look quite similar. The spectral structure in the 
spectrum of NGC 44 18 is, however, the product of strong ice 
and silicate absorptions (Spoon et al. 2001 ), whereas the spec- 
trum of M 83 is dominated by the commonly observed emission 
bands of Polycyclic Aromatic Hydrocarbons (PAHs). The spec- 
trum of Arp 220 shows characteristics of both: PAH emission is 
readily detected at 6.2 fim and 11.2 fim (and in the groundbased 
3 fim spectrum; Imanishi & Dudley 2000), while absorptions 
from water ice and silicates are found at 6.0 fim and 9.7 fim, 
respectively (Spoon et al. 2002). The strongest feature in the 
spectrum, peaking at 7.7 fim, has been previously attributed ex- 
clusively to PAH emission. On closer inspection, however, the 
feature seems to be a blend of a 7.7 fim PAH emission feature 
and a continuum peak resulting from absorption by ice on the 
short wavelength side and silicates on the long wavelength side. 

Close comparison of the 7.7 fim feature in the spectrum of 
Arp 220 with the 7.7 fim PAH feature in the spectrum of the 
starburst galaxy M83 (Fig.|3J shows that the 7.7 fim feature 
in Arp 220 is quite broad. Moreover, the ratio of the 7.7 fim 
peak to the 6.2 fim and 11. 2 /mi PAH bands is very large in 
Arp 220 compared to other sources (cf. Fig.[3J. The profile of 
7.7 fim PAH emission bands is, however, known to be very sta- 
ble over a wide range of integrated galaxy spectra (Helou et 
al. l2"0"00h . Peeters et al. (2002) found that Galactic ISM spectra 
show a similarly stable 7.7 ftm PAH feature. Only for evolved 



Fig. 2. A comparison of the ISO-PHT-S spectra of M83, 
Arp 220 and NGC 4418. While the spectrum of M83 is dom- 
inated by PAH emission bands, the spectrum of NGC 44 18 is 
dominated by absorption bands of ices and silicates. The spec- 
trum of Arp 220 shows characteristics of both. The spectra of 
M 83 and NGC 4418 have been scaled and offset. 



stars and isolated Herbig AeBe stars are the width and central 
wavelength of the 7.7 fim PAH feature known to vary (Peeters 
et al. I2002i . These stars are however very unlikely to domi- 
nate the spectrum of Arp 220. Therefore, there is no reason to 
assume that PAH features in Arp 220, if present, should have 
an intrinsically different shape than in any other galaxy spec- 
trum. Likewise, while the relative strengths of the PAH features 
are known to vary from source to source, no Galactic or extra- 
galactic sources show such extreme ratios as Arp 220 (Peeters 
et al. 2002). Alternatively, the large width of the 7.7 fim fea- 
ture might also be the result of a strong velocity dispersion 
within the PAH emitting gas (Av~5000 km/s). A similar ve- 
locity broadening is, however, not observed for the 6.2 and 
11.2/mi PAH features. Strong extinction as a cause for the 
unusual strength and width of the 7.7 fim peak can also be 
ruled out. This is best illustrated by fitting a starburst spec- 
trum (here: M 83) to the peak of the 7.7 fim feature and apply- 
ing foreground extinction (Fig. [3}- The weakness of the 6.2 fim 
PAH feature is reproduced well for a foreground extinction 
A(V)=105, but the spectrum beyond 8.5 yum and the 8.6/im 
and 11.3 fim PAH features are not. While the relative feature 
strengths depend on the adopted extinction curve (Sect. 3.2), 
no plausible extinction will be able to widen the 7.7 fim fea- 
ture. We here propose instead the weakness of the 6.2 fim PAH 
feature to be indicative of an unusually small contribution of 
the family of PAH emission features to the mid-IR spectrum of 
Arp 220 and the unusual strength and width of the 7.7 fim fea- 
ture to be the result of an underlying mid-IR continuum, peak- 
ing strongly near 7.7 fim. Evidence in support for this model is 
presented below. 
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Fig. 3. Comparison of the spectra of Arp 220 (thick grey line) 
and the central region of starburst galaxy M83 [black line). 
The M 83 spectrum shows strong PAH emission features at 6.2, 
7.7, 8.6 and 11.2 /an. In contrast, in Arp 220 the 6.2, 8.6 and 
11.2 /mi features are weak. The black dashed line illustrates the 
effect of strong extinction (A(V)=105; A(7.7/mi)=1.2) on the 
M 83 spectrum. Both versions of the M 83 spectrum have been 
scaled to match the 7.7 /mi peak in the spectrum of Arp 220. 



2.2. Broad 7.7 nm feature similar to Mon R2:IRS 1+2 

The shape of the 7.7 /mi feature in Arp 220 appears to be unique 
among a sample of more than 250 galaxies observed spec- 
troscopically in the 6-12/mi range (Spoon et al. in prep.). 
The feature has, however, a Galactic counterpart: the ISO- 
SWS spectrum of the combined line of sight to the sources 
MonR2:IRS 1+2 (Fig.0. MonR2 is a massive Galactic star 
formation region. Infrared observations show several compact 
sources and extended emission in the central region of the giant 
molecular cloud (Beckwith et al. I1976i . An elliptical ring en- 
closes two IR sources, IRS 1 and IRS 2. IRS 1, with a presumed 



spectral type of B0 (Downesetal. 1975; Howard et al. 1994J), is 
the exciting source of the compact H II region enclosed by the 
IR ring (Massi, Felli, & Simon [T985l . IRS 2 is still deeply em- 
bedded in the molecular cloud and most probably at an earlier 
stage of formation. Given this confused line of sight, the broad 
7.7 /mi feature in the spectrum of MonR2:IRS 1+2 may well 
be the result of the superposition of a strongly absorbed contin- 
uum, peaking at ~8 //m, and a 'normal' 7.7 /mi PAH emission 
feature. Fig.|4]shows the result of a crude decomposition of the 
MonR2:IRS 1+2 spectrum into the ISO-SWS spectra of the 
embedded protostar NGC 7538TRS 9 and the reflection nebula 
NGC 7023. The fit is quite good given the fact that the columns 
of ices and silicates vary a lot from one embedded protostar to 
the other. Note that the 8.6 /mi PAH feature is not suppressed 
by extinction, but instead stands out as a shoulder on the flank 
of the 7.7 /mi feature, filling in the blue wing of the 9.7 /mi sil- 
icate absorption feature. 
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Fig. 4. The line of sight to the embedded protostar 
MonR2:IRS 1 passes through the ultra compact Hn region 
MonR2:IRS2. The ISO-SWS spectrum of MonR2:IRS 1+2 
(black) hence shows features of both type of environments: 
PAH emission from the ultra compact Hil region and a 
dust and ice absorbed continuum from the embedded pro- 
tostar. Here we show a simple 2-component fit to the 
MonR2:IRS 1+2 spectrum, using the ISO-SWS spectra of the 
protostar NGC7538:IRS9 (light grey area) and the reflection 
nebula NGC 7023 (dark grey area) as fitting templates. 



2.3. Mid-IR continuum similar to NGC 441 8 

In Fig. [5] we compare the mid-IR spectra of Arp 220 and 
NGC 4418. The spectrum of NGC 4418 is dominated by strong 
silicate absorption bands at 9.7 /mi and 18 /mi (Roche et al. 
1986) and ice absorption bands at 6.0/mi (H2O), 6.85//m 
& 7.25 /mi (Hydrogenated Amorphous Carbons; HAC) and 
7.67 /mi (CH4) (Spoon et al. 2001 1. No mid-IR emission fea- 
tures, including the commonly detected 12.8 1 /mi [Neil] fine 



structure line, have been detected so far (Spoon et al. 2001; 
R. Siebenmorgen, unpublished TIMMI2 spectra). A first order 
estimate for the mid-IR local continuum of NGC 44 18 is ob- 
tained by fitting a power law through two feature-free "pure 
continuum" pivots at 8.0 /mi and 25 /mi (the dashed continuum 
in Fig.|5|i. Another, more conservative choice of continuum, as- 
suming "pure continuum" at 6.7 //m and 13 //m too, is shown as 
a dotted line in Fig. [5] 

The mid-IR spectrum of Arp 220 (Fig. [3 looks very sim- 
ilar to the pure absorption spectrum of NGC 4418, except for 
the presence of weak emission features due to PAHs (6.2 /mi, 
7.7/mi and 11. 2 / mi), 6. 91 /mi H 2 0-0 S(5) and 12.81/zm 
[Ne 11] (Sturm et al. 1 19961 Genzel et al. 1998 1. Using the same 
method as for NGC 4418, we estimate the local continuum in 
Arp 220 by fitting a power law through two feature-free "pure 
continuum" pivots at 5.5 //m and 25 //m (the dashed contin- 
uum in Fig.|5Jl- A more conservative estimate for the local con- 
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Fig. 5. Comparison of the mid-IR spectra of Arp 220 (multiplied by 90) and NGC4418. For Arp 220 the following spectra are 
plotted: in black the 5-16 //m ISO-CAM-CVF spectrum, in light grey the 5.6-1 1.4//m ISO-PHT-S spectrum and in dark grey 
the 17-22 yum UCL spectrum (Smith et al. 1989 1. ISO-SWS continuum points are marked by open circles, ISO-PHT photometry 
by filled circles and the IRAS 25yum flux by across. For NGC 4418 the spectra plotted are: in black the 5.6-11. 4 /mi ISO-PHT-S 
spectrum and in light grey the 8-13/mi and 17-23 yum UCL spectra of Roche et al. ( 1986 1. IRTF photometry (Wynn- Williams & 
Becklin 1993 1 is marked by filled circles, ISO-CAM (LW3) photometry by a triangle and the IRAS 25/mi flux by a cross. The 
positions of absorption and emission bands are indicated, as are several choices for the local continuum for each object. 



tinuum is obtained by including "pure continuum" pivots at 
6.7 yt/m and 14-15/mi as well. This results in the dotted con- 
tinuum for Arp 220. The dash-dotted continuum in Fig.|5]is a 
compromise between the two. Note that besides the obvious 
silicate absorption features at 9.7 /im and 18/im all three con- 
tinua imply the presence of an absorption feature due to water 
ice, which runs from 5.7 fim to 8.0 fim (Spoon et al. 2002 1. The 
emission features within this range, 6.2 /im PAH, 6.91 yum H2 
0-0 S(5), 7.7 yum PAH and likely (though not targeted by ISO- 
SWS) 7.0 yum [At ii], fill up the absorption partially or even turn 
it into emission. 

Given the close similarity of the mid-IR continua of 
Arp 220 and NGC 44 18 and the presence of strong silicate and 
ice absorption features in both spectra, the mid-IR spectrum of 
Arp 220 seems to be the superposition of a strongly absorbed 
continuum and a typical PAH-dominated spectrum. In the next 
Section we will test this hypothesis by decomposing the mid- 
IR spectrum of Arp 220 into a PAH-dominated spectrum and 
an absorbed continuum. 

3. Mid-IR spectral decomposition 

In order to test the superposition hypothesis, we have fit sev- 
eral combinations of a mid-IR absorbed continuum source and 
a PAH template to the observed Arp 220 mid-IR spectrum. Our 



method differs from existing methods (e.g. Tran et al. 2001) by 
using observed templates instead of model mid-IR continua. 
The latter usually do not take into account the complex radia- 
tive transfer effects due to ices and silicates that give rise to the 
exotic observed spectral shapes of the continuum sources and 
hence may fail to reproduce their spectra properly. 

3.1. PAH and continuum templates 

The absorbed continuum sources span a range of spectral 
shapes (Fig.|6^), with 6-12yum peak flux wavelengths ranging 
between 7.7 /im for 103344-2103 and 8.3 yum for the Galactic 
Center (SgrA*). All four sources exhibit a strong 9.7 yum sil- 
icate feature and three of them also show clear signs of ice 
absorption features. The silicate optical depths range from 
r sa ~ L9]for SgrA* and r sil >1.9forIRAS 00183-71 11 (Spoon 
et al. l2002l to T si] >3.5 for IRAS 03344- 2103 a nd r sil >3.7 for 
NGC 441 8. For the latter, Roche et al. (1986) claim an even 
higher value, r s ;i ~7. Judging from the steepness of the blue 
wing of its silicate feature (Fig.|6^), NGC 4418 may well have 
the largest dust column of all four absorbed continuum tem- 
plates. 

1 this is the aparent optical depth, measured from the ISO-SWS 
spectrum (Fig.|6^). Corrected for silicate emission along the line of 
sight the true value is 3.6±0.3 (Roche & Aifken ll9851 . 
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Fig. 6. PAH and continuum templates used in the decomposition of the observed mid-IR spectrum of Arp 220. Left panel: Four 
mid-IR absorbed continuum spectra. The spectra are shown at their instrumental resolution: R~90 for IRAS 03344-02103 and 
NGC4418; R~40 for IRAS 00183-71 1 1; R~l 000 for Sgr A*. The spectra have been normalized to their 7.5-8.5 fim peak fluxes 
and are sorted according to the central wavelength of their 7. 5-8. 5 fim peak flux. Right panel: Three mid-IR PAH-dominated 
spectra, shown at a spectral resolution R~90. The spectra have been normalized to their 7.1 fim peak fluxes and are sorted in 
order of increasing 9-12yum continuum. 



As PAH templates we selected the reflection nebula 
NGC7023, the Orion Bar star forming region and the cen- 
tral region of the starburst galaxy M83. The latter has been 
preferred over other nearby starburst templates, like M 82 or 
NGC253, because its 5-16 fim spectrum appears less affected 
by extinction than the other two galaxy spectra. Together, the 
three selected PAH templates (Fig.lU) are meant to cover the 
full range of spectral shapes from quiescent PDRs to intense 
star forming regions. Note the large difference in strength of 
the 9-12//m continuum between the spectra of NGC 7023 and 
the Orion Bar. The difference is attributed to the presence of 
hot dust in the Orion H II region, which is absent in the cold 
environment of a reflection nebula. In contrast, the PAH spec- 
tra are quite similar, except for a somewhat weaker 6.2 fim PAH 
feature relative to the 7.7 fim PAH feature in M 83 and a some- 
what stronger 11.2 fim PAH feature relative to the 7.7 fim PAH 
feature in the Orion Bar spectrum. 

3.2. The extinction law at mid-infrared wavelengths 

Unlike the extinction law at optical or near-IR wavelengths, 
surprisingly little is known about the shape and the general ap- 
plicability of the extinction law at mid-IR wavelengths. Results 
obtained for different lines of sight vary considerably. A good 
definition of the shape of the mid-IR extinction law is, however, 
important for a good decomposition of the Arp 220 spectrum. 

The differences among extinction laws appear largest at 3- 
8/im. Assuming a standard graphite-silicate mixture, a A~ lJ5 
power law fall-off is expected for this wavelength range (Fig.0 
Draine IT9891 Martin & Whittet IT9901 ISO measurements of 
molecular hydrogen towards the Orion Peak-1 as well as Hi 
recombination lines of compact H II regions support this model 



(Bertoldi et al. Martin-Hernandez et al. EES- H I re- 

combination line observations, probing the dusty, complex line 
of sight to the Galactic center (SgrA*), indicate, however, a 
nearly flat extinction curve between 3 and 8 fim (Fig.0 Lutz 
1999 1, with A,i/Ay a factor ~4 higher at 7 fim compared to the 
extinction law of Draine ( 1989). This would point to the pres- 
ence of larger, fluffier grains in the line of sight to the Galactic 
center than towards other Galactic sources (e.g. Kriigel & 
Siebenmorgen 1994 1. 



Beyond 8 fim, the 9.7 and 18 pm bands of amorphous sili- 
cates dominate the extinction curve. Both the shape (FWHM) 
and the strength of the bands (both Agj/Ay and Agj/Aig) are 
reported to vary between different lines of sight (e.g. Fig.[6k) 
and from observer to observer (Draine 1989). Here we will 
adopt the astronomical silicate profile of Weingartner & Draine 
(200U). 



In order to explore the effect of differences between mid-IR 
extinction curves on the spectral decomposition, we here define 
two extinction curves which should be representative for the 
range of extinction properties between different lines of sight. 
The extinction law, which we will refer to as Draine_local, 
combines the /T L75 power law fall-off for A <8.14//m with a 
silicate feature of strength A9 (,/Ay=0.06 typical for the local 
solar neighbourhood (Roche & Aitken [T984> . In contrast, the 
Lutz_gc extinction law is 'flat' from 3-8 fim and has a silicate 
feature strength of A 9 6 /A 1 /=0.14 (Lutz 1 1999k Both extinction 
laws are shown in Fig.0 



H.W.W. Spoon et al.: Mid-infrared spectral evidence for a luminous dust enshrouded source in Arp 220 



7 



t 1 1 1 1 — i 1 r 




3 5 10 20 

Wavelength [//m] 



Fig. 7. Comparison of two mid-IR extinction laws, overlayed 
on the PAH emission spectrum of reflection nebula NGC 7023 
(grey area). The Galactic center extinction law of Lutz ( 1999) 
is shown as a black line, the extinction law of Draine ( 1989) 
for the local solar neigbourhood as a dotted line. The extinc- 
tion laws are refered to as Lutz_gc and draine_local, re- 
spectively. 

3.3. Decomposition method 

For each continuum and PAH template combination we have 
explored a four parameter space for the best fit to the observed 
Arp 220 ISO spectrum. The four parameters explored are: 

- the contribution of the absorbed-continuum spectrum to the 
total spectrum 

- the contribution of the PAH spectrum to the total spectrum 

- the adopted extinction law: either draine_local or 
Lutz_gc 

- the amount of reddening (Ay) on the PAH spectrum 

For a 2-component model spectrum to give a good fit to the 
Arp 220 spectrum, the model spectrum has to reproduce several 
key spectral features which make the Arp 220 spectrum stand 
apart from other mid-IR galaxy spectra. These key features are, 
in order of decreasing importance: 

- the width of the 7.7 /mi feature 

- the absence of the 8.6/mi PAH feature in the red flank of 
the 7 .7 /mi feature 

- the strength of the 11. 2 //m PAH feature relative to the 
7.7 /mi feature 

- the strength of the 6.2 /mi PAH feature relative to the 7.7 /mi 
feature 

- the depth of the 9.7 /mi silicate feature 

- the profile of the 6.0/mi water ice feature 



3.4. Decomposition results 

In Fig.[8]we present the best fits to the 6-12/mi Arp 220 spec- 
trum. The fits clearly show that an appreciable contribution 
from a strongly absorbed continuum source (e.g. NGC 4418, 
or a stronger absorbed version of this spectrum) is required for 



a good fit to the observed spectral features. Compared to the 
extinction on the continuum source (r s n >3.7), the extinction 
on the PAH component is minor, ranging from T s ji=0.2 to 1.4). 
At these small obscurations, the choice of extinction law does 
not play an important role. The best fits further reveal that a 
direct measurement of the silicate optical depth from the ISO 
spectrum will be severely hampered by the presence of filled-in 
emission from the PAH component (see Fig.[8j, resulting in a 
serious underestimation of the true silicate optical depth. 

Other combinations of our continuum and PAH tem- 
plates fail to reproduce key features of the Arp 220 spectrum. 
Fits involving the continuum templates IRAS 03344-2103 and 
IRAS 00183-71 11, for example, fail to fit the width of the 
7.7 /mi feature, because their continua peak at too short a wave- 
length (Fig.|6^). The continuum of Sgr A*, on the other hand, 
does peak at the right wavelength, but its silicate feature is too 
shallow and too narrow to leave much room for any contribu- 
tion from the PAH component to the 8.5-1 1.5 /mi fit. Strong 
extinction will have to be imposed on the PAH component to 
minimize its contribution to the 8.5-1 1 .5 /mi fit. This, however, 
also minimizes the flux in the 11.2/mi PAH feature to beyond 
what is still consistent with the observations. A similar prob- 
lem exists for the Orion Bar PAH template. Its 10/mi contin- 
uum is the strongest among the three PAH templates (Fig.[6ji) 
and hence requires an appreciable extinction on its contribu- 
tion to any Arp 220 fit; too much for a good fit to the 11.2 //m 
PAH feature. PAH template NGC 7023 suffers from the oppo- 
site problem. Its continuum is too weak to produce good fits to 
the Arp 220 spectrum. 

Absorptions by ices play an important role in distorting 
the spectral shape of strongly dust enshrouded sources like 
NGC 4418 (Spoon et al. l2Wl Fig.^). Their impact is, how- 
ever, small for moderately absorbed spectra like the line of 
sight to Sgr A* (Fig.[U). On the other hand, as the abundance 
of ices is variable and changes from one Galactic molecular 
cloud to another, it is useful to assess the impact of an increased 
water ice abundance on our fits. We therefore ran tests with a 
Lutz_gc extinction law with five times higher water ice abun- 
dance and compared the best fitting parameters to those for the 
unmodified Lutz_gc extinction law. Only for those template 
combinations requiring high foreground extinction on the PAH 
template were the individual parameters found to change no- 
ticeably (but <10%). For all other combinations, including our 
best fits (Fig.[8}, the effects turn out to be negligible. Further 
experiments with an extinction curve supplemented with both 
the 6/mi ice and 6.8/mi HAC absorption features show that 
under these conditions successful fits can be obtained also with 
an (additionally obscured) Sgr A* continuum. This stresses the 
presence of ice absorptions in the continuum as a key require- 
ment for a successful fit, and one of the reasons why NGC 4418 
produces the best results in our original fits. 

4. Discussion 

Our successful decomposition of the 6-12/mi ISO-PHT-S 
spectrum into a strongly absorbed continuum and a weakly ab- 
sorbed PAH component confirms our initial suspicion that a) 
its exotic mid-IR spectrum resembles a blend of the spectra of 
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Fig. 8. Four 2-component fits to the Arp 220 mid-IR ISO spectra (ISO-PHT-S: dark grey histogram; ISO-CAM-CVF: dark 
grey line). In each panel the dark grey area represents the contribution of the continuum component and the light grey area the 
contribution of the PAH component. The PAH component is the same in all panels: the spectrum of the starburst galaxy M83. 
The absorbed continuum component differs between left and right panels. In the left panels it is NGC4418 as observed; in the 
right panels it is NGC4418 as seen through an additional screen of A(9.6/im)=2.5, equivalent to a factor ten more attenuation 
at 9.6/im. The fractional contribution of each fit component to the total 6-12//m fit is stated in brackets behind the name of the 
template. Also indicated is the amount of 9.6 fim screen extinction on the PAH component. 



NGC 4418 and M 83 (Fig.|3 and b) that the only difference be- 
tween the mid-IR spectra of Arp 220 and a strongly ice- and 
dust-absorbed source like NGC 4418 is the additional presence 
of PAH emission features in the spectrum of Arp 220 (Fig.[3}. 

4.1. Identification of the spectral components 

The large difference in obscuration of the two spectral com- 
ponents of our fit indicates that these most likely represent 
two spatially separate components. Recent high angular res- 
olution N-band spectroscopy of the nuclear region (Soifer et al. 
2002 1 show the 11.2/mi PAH emission and the ll-12yum con- 
tinuum emission to have clearly different distributions. While 
the absorbed-continuum emission clearly peaks on the two nu- 
clei, the PAH emission extends over a far wider area and does 
not peak on either nucleus. We therefore associate our absorbed 
continuum component with the two nuclei and our PAH com- 
ponent with the region in between and around the two nuclei. 

We estimate the infrared luminosity associated with the dif- 
fusely distributed PAH component from the observed 6.2 /jm 



PAH emission feature by assuming a generic conversion ra- 
tio between L(6.2yum PAH) and L(IR) and taking into account 
the weak obscuration on the 6.2//m PAH feature as indicated 
by our best 2-component model fits. We derive the L(6.2jt/m 
PAH)/L(IR) ratio from our sample of ~70 mid-IR ISO spec- 
tra of normal and starburst galaxies and obtain a value of 
0.003+0.001. Assuming this conversion factor applies also to 
the conditions in the ULIRG Arp 220 and applying it to our 
four best fitting template combinations (Fig.[S} we find the in- 
frared luminosity associated with the diffuse PAH component 
to amount to 1 .2-2. 1 10 1 1 L ; 9-15% of the total infrared lumi- 
nosity of the system. A comparable result is obtained from the 
peak flux density of the 7.7 yum PAH feature, using the empiri- 
cal conversion factor S(7.7 p.m PAH)/F(IR)=1Q- 1L84 Jy W" 1 m 2 
for starburst galaxies, determined by Lutz et al. ( 2003 1. Taking 
into account the weak obscuration on the 7.7 yum PAH fea- 
ture (ranging from A(7.7yum)=0.07 to 0.49), the implied in- 
frared luminosity is 0.7-1.8 10 11 L ; 5-13% of the total in- 
frared luminosity of the system. The results from both meth- 
ods are in complete agreement with the value derived by Soifer 
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et al. ( 2002 ), who infered a PAH-associated infrared luminos- 
ity of 1.2 10 11 L Q (9% of the total infrared luminosity) from 
their N-band spectrum. Based on the small scatter between the 
three different methods, we conclude that the infrared luminos- 
ity associated with the diffusely distributed PAH component in 
Arp 220 amounts to 5-15% of the bolometric luminosity of the 
system, with some uncertainty due to the trend towards a larger 
FIR/PAH emission ratio with increasing interstellar radiation 
field intensity in galaxies (e.g. Dale et al. 2001 ). 

The literature reports the detection in the nuclear region of 
several near- and mid-IR recombination and fine structure lines 
from star formation. While no spatial information is available 
as to the origin of the mid-IR lines (Sturm et al. [1996 1, both 
the near-IR recombination lines (Armus et al. [1995 Larkin et 
al. Q995 ) and the radio recombination lines (Anantharamaiah 
et al. 2000 1 are mostly concentrated towards the two nuclei, 
rather than following the diffuse PAH component. Extinction 
estimates range from A(V)~10 in the near-IR (Armus et al. 
[T9951 Larkin et al .fT995l up to A(V)~40-60 in the mid-IR and 
radio (Sturm et al. 1996 1 Anantharamaiah et al. 2000 1, with not- 
icable uncertainties e.g. due to the measured line fluxes. This 
suggests that the starburst activity traced by the emission lines 
includes an intermediate obscuration component that is more 
concentrated on the nuclei than the PAH emission. Two factors 
suggest this component includes intermediate obscuration parts 
of the regions around the two nuclei rather than the deepest em- 
bedded parts: A(V)<50 appears still moderate compared to the 
extremely obscured NGC 4418-like continuum, and the energy 
budget of this component stays a factor ~2 below the bolomet- 
ric luminosity of Arp 220 (Genzel et al. 1998 1. The infrared 
emission lines may thus give a partial view of the circumnu- 
clear absorbed region, but very likely not a complete one. 

4.2. The nature of the nuclear power sources 

Given the relatively small fraction of the bolometric luminos- 
ity associated with the diffuse PAH emission (5-15%), the bulk 
of the luminosity of the Arp 220 system must be associated 
with the absorbed-continuum component from the two nuclei. 
To be responsible for only ~40% of the 6-12/im luminosity 
but for 85-95% of the total Arp 220 infrared luminosity, these 
nuclei must, hence, be deeply enshrouded indeed. Given the 
large amounts of molecular material detected in the nuclear en- 
vironment (~1O 1() M half of which is in the disk enclosing the 
two nuclei; Scoville et al. 119971 Sakamoto et al. [H>99l, this 
comes as no surprise. Unfortunately, strong obscuration erases 
source-specific spectral signatures, not only at X-ray, UV, op- 
tical, near- and mid-IR wavelengths, but, depending on the ab- 
sorbing column, also out to far-IR and, possibly, sub-mm wave- 
lengths. 

Based on the few identified spectral signatures at hand, we 
consider two possibilities for the nature of the nuclear power 
sources. First, each of the nuclei may contain a deeply embed- 
ded, extremely dense and luminous stellar cluster - a super-star 
cluster containing some 10 6 massive O stars within a region 
less than lOOpc in size. Less extreme deeply embedded clus- 
ters have been discovered in the starburst galaxies NGC 5253 



and He 2-10 (Gorjian et al. EUoTI Vacca et a!. 12UU2l A stel- 
lar nature for the nuclear power sources in Arp 220 is fur- 
ther supported by the 18 cm VLBI observations of Smith et 
al. J1998lfT9^9T >. which show a dozen or so sources scattered 
over the two nuclei. These "knots" are consistent with free- 
free emission from luminous radio supernovae expanding in the 
dense (circum)nuclear environment (Smith et al. 1998 1999 1. 
Second, we consider two deeply embedded AGNs. Despite the 
lack of AGN features in any waveband, a deeply embedded 
powerful AGN at the center of each nucleus cannot be ruled 
out. The column density required to block the AGNs from de- 
tection by BeppoSAX in hard X-rays is N H2 = 10 25 1 -10 25 3 cirT 2 
with a large covering factor (Iwasawa et al. 12001b . The non- 
detection of a high-T/, radio core in the 1 8 cm VLBI maps of 
Smith et al. (1998l ll999l would then imply the AGNs to be 
radio-quiet or strongly free-free absorbed. 

4.3. Modeling of the nuclear continuum 

In order to test whether the low mid-to-far-IR continuum ratio 
in Arp 220 may be attributed to strong dust obscuration on the 
emission of two identical, deeply buried, energetically domi- 
nant, nuclear sources, we used the dust radiative transfer code 
of Siebenmorgen et al. ( fT999l l200lt to model the 2-3000 pm 
nuclear spectrum. 

In the model we assume a cluster of OB stars, with stellar 
densities similar to those inferred for ultra-dense H 11 regions 
(UD H 11 regions; Vacca et al. l2002> . to reside at the center of 
each nucleus. For a given luminosity of 6 10" L , the stars will 
occupy a sphere with radius r=25 pc. Adopting a nuclear gas 
mass of 1O 9 M , spherical symmetry, constant gas density and 
a silicate optical depth of T s ji=4.6 (consistent with the silicate 
optical depth for our best fitting absorbed continuum source 
NGC 4418), we are forced to place the outer radius of the dust 
shell as far out as ~325 pc in order to lower the dust tempera- 
ture to values consistent with the observed 12-50/mi SED. The 
FWHM of the 1 1 pm light profile is ~20 pc, consistent with the 
observations of Soifer et al. (2002i. As the temperature in the 
dust shell does not drop below 40 K at the outer edge, we con- 
sider an additional cold dust component to account for the cold- 
est dust in the system. This component likely accounts for the 
far-IR/submm emission associated with the diffuse PAH com- 
ponent and for far-IR/submm emission from the nuclear region 
not covered by our 'simple' model. Here, we characterize the 
cold dust component by a grey body spectrum of 30 K with a 
dust emissivity index of 1.8 and a dust optical depth of 1 at 
100 pun. Both the model spectrum of the nuclei and the grey 
body spectrum are shown in Fig.[5] The Figure also shows the 
assumed contribution of the diffuse PAH component (M 83; see 
Sect. 3.4) to the Arp 220 mid-IR spectrum. 

We conclude from our modeling that the low ratios 
S(6^m)/S(60yum) and S(6yum)/S(100jt/m) can, in principle, 
be explained by the effects of strong dust obscuration on 
two deeply buried, energetically dominant, nuclear sources. 
However, this model predicts a hydrogen column of only 
~10 23 2 crrT 2 and hence favours stellar heating, as this column 
is a factor ~100 too small to account for the upper limits on the 
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Fig. 9. The 2-3000 fim spectrum of Arp 220 (grey) overplotted with a three component fit. The black dashed curve represents the 
model SED for the two identical nuclei. The black dotted curve represents grey body emission from 30 K cold dust. The black 
dash-dotted curve is the sum of these two components. The black continuous curve is the sum of the model SED and the weakly 
absorbed diffuse PAH component (M 83) from our 2-component decomposition (Sect. 3). 



hard X-ray flux if the total luminosity were dominated by AGN 
activity (Iwasawa et al. 2001 1. 

Note that if any deeply hidden AGN (N H >10 25 cm -2 ) were 
present, the covering factor of the obscuration would have to be 
large in order to be consistent with the lack of reflected X-ray 
light measured by BeppoSAX (Iwasawa et al. 2001 ). The AGN 
would hence strongly contribute to the far-IR continuum, but 
could not be responsible for the mid-IR continuum. 

4.4. Mid-to-far-infrared spectral characteristics 

Arp 220 is an outlyer in many spectroscopic diagnostic dia- 
grams linking mid-IR to far-IR quantities. In all cases Arp 220 
stands out by having a low ratio of the mid-IR characteristic 
with respect to the far-IR characteristic: 

- Already in the Bright Galaxy Sample of Sanders et al. 
dl988> . Arp 220 stands out by the smallest S 12/S60 ratio. 

- Similarly, Arp 220 has the lowest S5.9/S60 ratio among the 
larger sample studied spectroscopically in the mid-IR by 
Lutz et al. ( fT998b . 

- The ratio of F(6.2yum PAH)/F(FIR) is lower than for any 
other galaxy in our database of more than 250 ISO mid-IR 
galaxy spectra. 

- The ratio F([Nen])/F(FIR) is a factor three lower than for 
the average starburst galaxy (Genzel et al. 1998 1. 

- The ratio F(7.7yum PAH)/F(850/vm) is the lowest among a 
sample of normal and ultra-luminous galaxies (Haas et al. 

i2oon> . 



Our spectral decomposition offers a simple explanation for 
the above observations. A typical starburst spectrum associ- 
ated with the extended PAH component will contribute strongly 
to the mid-IR but weakly to the far-IR, while the domi- 
nant deeply enshrouded nuclei contribute weakly to the mid- 
IR and strongly to the far-IR. The resulting combined 3- 
1000 /mi SED is hence mostly starburst-like in the mid-IR 
and dominated by cold dust emission from the enshrouded 
nuclei in the far-IR. Ratios of mid-IR to far-IR quanti- 
ties, like F([Nen])/F(FIR), F(6.2/im PAH)/F(FIR) or F(7.7 pm 
PAH)/F(850yum), will hence all be systematically lower than 
the values typically found for starburst galaxies. Interestingly, 
Luhman et al. (2003 1 advocate a similar scenario as a contrib- 
utor to the low F([C n])/F(FIR) ratio of Arp 220. 

Other galaxies may exist with even more strongly obscured 
nuclei, emitting an even smaller fraction of the nuclear 
luminosity in the mid-IR than Arp 220 does. These galaxies 
would hence look starburst-like in the mid-IR (contributed by a 
[weakly obscured] circumnuclear starburst) but would be char- 
acterized by a very strong cold dust continuum in the far-IR. 
Galaxies of this type, may be recognized spectroscopically by 
their low S5.9/S60, S5.9/S100, S5.9/S850 or F([C n])/F(FIR) 
ratios and their starburst-like 6.2//m-PAH line-to-continuum 
ratios. Perusal of our ISO spectral database has led to the 
identification of several candidate galaxies. Arp 220 may thus 
be a local and less extreme template of the class of SCUBA 
sources seen in the Hubble Deep Field. 

Arp 220 is not the only galaxy with spectral structure in the 
6-12//m range reminiscent of strongly modified PAH bands. 
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Other examples are Mrk231 and most of the galaxies in the 
sample of Tran et al. (©Oi l. Some of these spectra show, in 
addition to strong silicate absorption longward of 7.7 /mi, clear 
signs of water ice absorption shortward of 7.7 yum. Depending 
on the strength of the PAH emission features these galaxies 
have been classified as class I, II or III ice galaxies (Spoon et 
al. 2002 1. Like Arp 220, their spectra may well be the result of 
strong extinction on one or more spectral components. 

5. Conclusions 

We have shown that the 6-12yum spectrum of Arp220 is not that 
of a scaled-up typical starburst galaxy but contains a 'normal' 
starburst component characterized by PAH emission features 
plus a highly obscured dust continuum with ice and silicate ab- 
sorption. Attempts to decompose the spectrum using a variety 
of extragalactic and Galactic template spectra yields a best fit 
in which a typical starburst, represented by M83, contributes 
~60% and and an ice absorbed continuum galaxy, represented 
by NGC4418, -40% of the 6-12/mi luminosity. An impor- 
tant result in relation to previous studies is our conclusion that 
the pronounced emission feature peaking around 7.7 fim is a 
blend of PAH emission and a broader feature in the contin- 
uum caused by ice absorption at shorter and silicate absorption 
at longer wavelengths. We tentatively conclude that the PAH 
emitting component is only weakly absorbed and arises in the 
extended region imaged at higher resolution around 10//m by 
Soifer et al. (|2002 1 whereas the absorbed continuum is asso- 
ciated with one or both of the compact nuclei. This extended 
starburst component contributes only 5-15% of the total lumi- 
nosity with the bulk emitted by the heavily obscured nuclear 
component(s). One possibility is that this luminosity is gener- 
ated by starburst activity occuring in a higher density environ- 
ment than found in lower luminosity starburst galaxies due to 
the larger quantity of molecular gas and dust funnelled to the 
center by merging of the two nuclei. Due to the high extinction, 
it is not possible with the available data to exclude that AGN(s) 
also contribute some or all of this luminosity. Based on the up- 
per limits for hard X-ray emission (Iwasawa et al. 2001 ), how- 
ever, Arp 220 would need to be the most highly obscured AGN 
known. 
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